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ABSTRACT 


This thesis explores the possibilities of developing fast grid-based parallel algortthinis to 
solve the Weighted Region Least Cost Path problem. Two complementary steps have been 
mdertaken. First. an efficient sequential algorithin to solve the above problem was devel- 
oped. The algorithm is a modification of a Gauss-Seidel-like algorithin for obtaming the 
minimnunm costs. The most salient feature of the algorithm is the reduction of the number of 
nodes and edges in cheaper regions of the gud. The reported experimental results ascertam 
the superiority of this algorithm with regard to computer running time ata modest reduction 
inthe accuracy of the obtamned solution. Parallel implementations of grid-based algorithms 
were studied. A simple grid-based variant was fimplemented on a network of Transputers. 
The overall approach employed could be used to develop a parallel version of the above 
sequential algorithm on a Transputer network, combining both advantages of efficiency 


and parallelism. 
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I. INTRODUCTION 


A. GENERAL 

Numerical path planning has been studied quite a bit in 
the last few years. One problem of numerical path planning 
involves finding the optimal path from a given starting point 
to a goal point through a plane that has been subdivided into 
weighted regions. This problem is known as the Weighted Region 
Least Cost Path (WRLCP) problem. The best path can be 
minimizing some cost (e.g. energy or time). 

Several approaches can be used to tackle this WRLCP 
problem. Some of these approaches and their pros and cons are 
presented in chapter II. 

One of the most significant constraints has been the 
ability to evaluate in real time the optimal path through a 
weighted region. Traditional sequential algorithms can 
quickly become overloaded with the sheer number of 
calculations required for a realistic problem. Consequently, 
several algorithms for solving this type of problem have been 
researched in recent years. These algorithms have, to a 
varying degree, potential for parallelization. Parallel 
algorithms have a tremendous advantage in speed, but also may 
open the door to new problems, especially in mapping, 
scheduling and cloroneeinarsne, mare the She eng Shane parallel 
activities. Multicomputer networks are no exception. The 
principal goal in parallelizing an algorithm is to sustain a 


close-to-linear speedup in processing time. This is no trivial 


task, since the processor communication reduces the speedup. 
This, and the fact that parallelizing an algorithm may 
increase its overall processing time (summed over the total 
network of processors) make this problem challenging. 

We investigated in this thesis some of the basic approaches 
with their potential for parallelization, and implemented one 
of these algorithms, first sequentially and then in parallel 


on a Transputer network. 


B. THESIS OUTLINE 

Chapter II introduces an analysis of the different 
parallel path approaches which have been investigated 
recently. The algorithm which appears most promising for 
solving the weighted region least cost path problem, which is 
a grid algorithm, has been selected for further study. In 
chapter III, implementation and some modifications of this 
algorithm is presented with an illustrative example. Chapter 
IV provides a description of the Transputer network and ADA. 
In chapter V, this approach has been developed into a parallel 
algorithm and implemented on an INMOS Transputer network using 
ADA. The algorithms ability to use parallel computing to solve 
arbitrary WRLCP problems has been investigated. Finally, 
these results will be compared with the results of single 
processor, and with previous results on a multicomputer 
network. Conclusions and recommendations for further research 


are offered in chapter VI. 


It. BACKGROUND 


We overview in this chapter some of the basic approaches 
for solving the WRLCP problem and we investigate their 
potential for parallelization. 

A. SINGLE SOURCE SHORTEST PATHS (SSSP) 
1. Dijkstra’s Algorithm 
The algorithm for SSSP presented in Figure 1, which was 
@evelooed by Dijkstra im 1959, as the startinggpoint for this 
problem [MAN89]. 


Algorithm Single_Source_Shortest_Pathts (G,v); 
a Gi: G=(V, E) (a weighted directed graph) , and v (the source vertex). 
Output for each vertex w, w.SP is the length of the shortest path from v to w. 
(all lengths are assumed to be nonnegative. } 


eeu 

or all vertices w do 
w.mark := false; 
w.SP: =infinite: 

end loop; 


veoie = 0; 
while there exists an unmarked vertex do ——si*=n 
let w a an pon aree vertex such that w.SP is minimal; 
w.mark := tru 

for all a2 (w,z) such that z e pningiced do 

if w a ee ae en 
+ length(w,z) 
PCa 
end loop; 


oh ay 





Figure 1: Dijkstra’s Algorithm 


In Figure 2, a small example 1s presented to 
demonstrate the algorithm. The first line includes only paths 
of one edge from v (the source). The shortest path is chosen, 


in this case, leading to vertex a. The second line shows the 


update of the paths including now all paths of length one from 
either v or a, and the shortest path now leads to c. A new 
vertex 1s added in each line, and the current known shortest 
paths from v are listed to every vertex. The underlined 
distances are those that are known to be the shortest. The 
algorithm keeps adding new vertices to the selected list until 


all vertices are added. 





Figure 2: An Example Of The SSSP Algorithm 


The algorithm can be easily extended from directed 
meomundrrected form. 
2. Parallelization Of SSSP 
Dijkstra’s algorithm has been parallelized by Moore 
[OUI87]. He devised two parallel algorithms. The first 
algorithm makes the for loop (in Figure 1) parallel, which 
explores the outgoing edges from a given vertex. the second 
method is to parallelize the while loop (in Figure 1); that 
is, at any one time in the execution of the algorithm there 
are probably many vertices in the queue. The parallelizability 
@iertne Eirst method 1S restricted by the number of edges 
outgoing from each vertex. On the other hand, the second 
method performs larger tasks, that produces good speedup. More 
detailed information can be found in the given reference. 
In the parallel algorithm based on the second method, 
a queue is used. That queue is initialized with the source 
point, and then a number of asynchronous processes are 
created. Each of these processes goes through the steps of 
deleting a node from the queue, examining its outgoing edges, 
and inserting into the queue the nodes to which shorter paths 
have been found. In Figure 3 an example is presented in which 
the number of nodes in the queue shows actually the number of 
processes that can be parallelized. Distances are kept until 
they are reached from another direction. When there are more 


than one edge reaching to a node, minimum cost is chosen. 





Figure 3: Parallelized SSSP Algorithm Example 
B. RAY TRACING 


1. Snell’s Law of Refraction 
Snell’s law [HEC87] defines the path a ray of light 
takes as it passes from one medium to another. The ray is 
refracted across the border between the media according to the 
following equation, 


n,sin8, = n,sin8, 


where nN; Mp are the indices of refraction and the angles are 


of incidence and refraction respectively (Figure 4). Fermat’s 
principle which implies Snell’s law states that light follows 
a path between two points such that it takes the minimum time. 
It can be proven (for example see [RIC87]), that a similar 
principle govern the path a particle takes across two regions, 


in which the speed of the particle is uniform in each region. 


Therefore, we can apply the principles of optics to solve the 
WRLCP problem by assuming regions as optical media and weights 


Be imecitees, G@ refraction. 


medium 1 
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Figure 4: Snell’s Law Of Refraction 


2. Implementation 

The main step of this implementation 1s to shoot a ray 
from the source and see where it 1S going to land. We keep 
felng Eta wter Intttal raysewath different angles until some 
of the rays hit the goal. In essence, the boundaries border 
the weighted regions, the index of refraction depends on the 
weights and the rays refract on the border to border until it 
reaches the goal. At every boundary of two regions the ray 
obeys Snell’s law of refraction. Among those rays that hit the 
Geaepelme we  Olralin the WRLCP. In Figure 5, at the source 
point the short rays are just to give an idea about different 
set of angles and continuing three rays are given to make a 


clear example. 





Figure 5: Ray Tracing 


Through different weighted regions, as the path Sklmnop 
misses the goal, the minimal paths SabcdeG and SfghijG hit it, 
and one of them is chosen as WRLCP after comparison. 

3. Parallelization 

The Snell’s law based algorithm can be easily 
parallelized by assigning a different set of angles to each 
processor. If a solution 15 found, 1€ 15 optimal. @hercoesure: 
finding the intersection point of a ray with a triangle is not 
large. However, the algorithm suffers several drawbacks: The 
ray inversion, possible presence of blind regions, and use of 
expensive trigonometric functions. More details can be found 


ine (REGS r ane 


C. EDGE SLICING 

1. Implementation 

In this technmaque, every edge is divided into a number 

of segments of equal length. Some points are placed 
equidistantly on every map edge in the triangulated plane. A 
graph 1S constructed by connecting every two points on two 
edges belonging in the same triangle. The distance between two 
consecutive points is made proportional to some function of 
the costs of the two regions separated by that edge. A graph 
1S constructed whose nodes are the original triangles’ 
Vertmees plus the points which divides the edges into 
segments. The edges of this graph are the original triangles’ 
edges plus the lines connecting every two non-colinear points 
in the same graph. In Figure 6, edge slicing is shown for the 
case where an edge is divided into two equal segments. 

2. Parallelization 

After the edges ina graph are sliced, the WRLCP problem 

is reduced to a SSSP problem, and the technique used for 
Baralielizing SSSP@@an be invoked. A variant of this algorithm 
1s shown in [KIN91]. 
D. GRID ALGORITHMS 

Another approach for approximating the WRLCP problem is to 
model the terrain as a grid. Simply, a grid is laid over the 
terrain. The map 1S divided into equidistant grid points. The 
Weights in the regions are tramsferred as edge costs. Figure 
7 1llustrates how a node is connected to different number of 
neighbors (we implemented a node with 4-neighbors, but it is 


easy to extend the implementation for more neighbors). 





Figure 6: Edge Slicing (2-Point) 


Finally, the shortest path analysis 1s performed on the graph. 
Two classes of grid algorithms could be used. These are 
explained below. 
1. The Wavefront Variant 

Starting with the source point at) time 92ero see 
progresses one step every time unit in all directions adding 
the appropriate nodes (these nodes which can be reached in the 
earliest at this time step) to a wavefront depending on the 
edge weights and direction of propagation. The wavefront keeps 
advancing every time unit until it hits the goal point. It is 
obvious that the progress in the inexpensive areas 1S bigger 


than expenSive ones, since the léngth of step ts preeporueteme 
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Figure 7: Neighbors In Grid Algorithm 
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to cost and time. Figure 8 shows an example in which the bigger 


Strides take place in the low cost areas. 
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Figure 8: Wavefront Technique 


One wavefront sweep will find the shortest path. But 
such shortest path may need a large number of time steps 
depending on edge weights. The number of nodes on the 
wavefront to be processed at every perme Iieime valle, some 
causes difficulty in scheduling “tmem "on "a  Geeoeurer-—« 
computer. 

2. Relaxation-Based Approach 
The reason why this approach is called “relaxation- 


based” 1s because it 18S similar to Gauss-Seide! 1§eraticn 


eZ 


class of computation in Partial Differential Equations 


[HAB87] . 


oee(m + 1) = Cee Ole) eG.) et) 


He, ey) Capea in + Lede, Cm) | 


Pe cme sp erian) Shows value "Or Grid polne in row 7, column j of 


step m. 


The bracketed term indicates the change that occurs 
eee cache ererataon as cyeyim) is updated to cy,4 (m+l). 


Similarly, in our implementation, as shown in Figure 9, in 
each iteration, a node is updated. The cost of node is 
minimized by comparing it with the neighbors according to the 
following pseudocode: 

Ceurranni7 2) =entn £ Coig(i,j3), ¢(1,j3+1) + north_weight 
c(itl,j) + east_weight 
c(i,j-1) + south _weight 

c(i-1,j) + west_weight ] 

We decided to implement the relaxation-based grid 

algorithm because of the straight-forward solutions on a 
Srserioulced COomouter - 
a. Implementation 
The grid graph is represented as a two-dimensional 
ae eOmmEeecercnm Data Structure andwinformmation carried by 
every node 1s shown in Figure 10 and the variables are 


described below. 
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(1,j+1) 


north_weight 


ct wes i eas | 


(3057) 5), (1,3) (eel) 


south_weigh 


is (1,j=1) 





Figure 9: A Node Described With Its Neighbors 
And Their Relations 


Current_cost: The beginning and updated cost value 
of the node. 

Old_cost: Old_cost is the previous iteration cost 
value of the node and used for determining the change in cost. 

N,E,S,W: North, East, South and West neighbors. 


Weight: The weight on edge toward that direction. 
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(1,n) (t+1,n) 


GRID_POINT (record) 


urrent COs 


Vid cost 


ditice— 
| disionce 
el 


& o ny 


3B Tai 


(Ey (eT) 
*Direction of minimum path 
passing through node 


GRID 


(2-d array of records) 





Figure 10: Data Structure 


Distance: The step between two nodes (here it is 
always equal to one, but it changes in the modified 
Seugorithm) . 

Direction: Semi Neotel. East”, South”, “West ”) 
and used for showing the path in the output. 

The pseudocode for the algorithm is presented in 
Figure 11. 

We obtained a data file which has 6400 raw data 
values representing the heights of 80*80 grid approximation 
of region. We constructed costs for the edges of the grids 


based on some function of these heights. For border nodes, 


i> 


Algerithm: To find the Weighted Region Least Cost Path 


Given: Two dimensional grid of points: GRID 
SOUTEGE Polme (1s ,75) 
Goal petphsa 4 woe jG. 
Threshola er 


Qutput;: COSt of minimization fErom source to gaa: sees 
The minimum path from source to goal: 


(1s,3s) ld , jee 
procedure: 


Initialize 
Read the data file 
for all nodes loop 
calculate weight (1,)) 
end loop 
xx --At this point code will be inserted in the 
~-modified version of the algorithm (see Fig.14) 
flag “= trite 
while flag loop 


flag = false 
for all nodes loop 


find minimum value 
change = node(i,j3) .old_cost-node(i,4) -currenemeea 
1f change > T then 
raise flag 
end if 
end loop 
end loop 


output the WRLCP 





Figure 11: The Psuedocode For Relaxation-Based Algorithm 


edges leading out of the grid where assigned very high 
positive values for weights, and minus one to distance. For 
all nodes we initialize the costs to a maximum positive number 
in order to use them in comparisons for finding the minimum. 

To output the Least Cost Path, we utilize a stack 
Since the path traces back from the goal to the source point. 


So the output is displayed from Source to geal@acmces 


ime 


direction and distance (distance is equal to one in the 
straight-forward algorithm) are calculated at each step. 
A small-sized (4*4) example is presented in the 
Bagure IZ £or clear understanding. 
b. Parallelization 
In the relaxation-based approach, solutions can be 
straightforward on a vector computer or a on 2-d mesh of 
processors. But also there are some unattractive features such 
Peee@e algoriram 1S intrinsically nonmoptimal (that is true for 
the wavefront technique, too). These algorithms do not reward 
the areas with the low cost by reducing the computations 
there. If the relaxation approach 1s used on a realistic 
Mimoe> Of Drecessors, It 1S not known how to partition the 
Sempuealions. FOr instamee, only one of the processors has the 
source point and starts computing, the computations of the 
rest of the processors for taking the minimum of infinite 
values are useless until they get the border values from that 


processor. 


Step 1. The heights are read from the file 


and loaded to nodes. 


Step 2. The weights on the edges are calculated 
using some functions 


(e.g. abs(1024 - 1023) + 1 





Figure 12.a: Example For Relaxation-Based Algorithm 
The Values Inside The Nodes Are The 
Heights 
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feoneimuca rrom Paoqure t2 7a) 


Step Sesounce™po2nt current ,cost = 


The others.current_cost = 





Figure 12.b: Example For Relaxation-Based Algorithm 


icmiimes= Im mie Nedes Are The Initial 
GastsuerFrom The Source. 


lle, 


(continued from Figure 12.b) 


See pm 


Figure 12.c: 


This 1s the first iteration. In the 
latter iterations continue. 

14 becomes 12 after first iteration 
Of COSE minimiZamaen- 

Thick arrows show the minimum least 


cost path 





Example For Relaxation-Based Algorithm 
The Values Inside The Nodes Are The 
Minimum Costs From The Source 
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IIt. A MODIFIED GRID ALGORITHM 


A. INTRODUCTION 


ireweconplexity Of the straight-forward grid algorithm in 


worst-case 1S Oa rEoumemrn x 7 Crla. tn the worstmease, the 
WRLCP passes through all the nodes. Since the algorithm does 
not take the advantage of the low cost regions, we modified 
mie ral gorithm to remedy this defficiency. 

B. VARIABLE GRID SIZE APPROACH 

We can use relaxation to obtain the shortest path from the 
Seurce tO different points without updating the nearest 
neihgbors as in the classical technique. We update the 
farthest neihgbors (two horizontal and two vertical in case 
of a 4-neihgbor technique) whose weighted distance is less 
than or equal to a prespecified constant. Determining such 
neighbors is done at the initialization phase. Rather than 
having the wavefront advance by a fixed distance, this 
algorithm has it advance by a fixed weight, thus bigger 
strides are made in cheaper regions. 

Information that decides which neighbor to propagate to 1s 
shown in Figure 13. It assumes that every grid point will have 
four records (N,E,S,W), and every record will have two 
components: the number of grid points to the neighbor to be 
updated (one in the classical method), and the weighted 


distance to the neighbor. 


IEA 


Algorithm: I.To make horizantal search and elimination 
TI.To decide which neghbor to propagate 


Given: Two dimensional grid of points: GRID 


(all edges 
Source point: (2s), 73s) calculated) 
Goal point = sever 1c) 
Stride 
Qutput: A searched grid and neighbor for node(i,j) to 
propogate 
procedure horizantal search and elimination 
for all rows 1 loop 
J = 
while j < jmax loop 
J = Jo 
loop 
z£ weight from (1S5S,37S) tomiig sq) soetben 
I) =e aie 
end if 
until (weignt frome(is, ys) tee ye — or 
or ({j = Jmaee) 
if weight from (l1sejs) to o(tg7g6) © => tuner 
Je= ac! 
end if 
delete all horizantal edges between (1,jg) and 
| lh) 
(1,J9)-east_weight = (1,jJ).west_weight 
= weight difference 
(1,J9).east_distance = (1,]).west_distance 
=o 
end loop 
end loop 


procedure: propogation 


for all nodes (1,)j) loop 
for north loop 
while {(1,]).north_weight + (1,]) .north_weight 
(1,3) .north_distance} < stride loop 


(1,]).north distance: =(19 49) = nore aistene oma 
(1,)]) .north_weight:= old value + 
(1,j3+(1, 3) .north_distance) .north_weight 
end loop 
end loop 
end loop 


repeat for other directions 


Figure 13: The Pseudocode for Modified Grid Algorithm 


ae, 


The procedure propagation assumes that the distance 
between successive grid points has been already computed; 
otherwise a more efficient strategy could be utilized. 

ii miecwelecd=tolaxarlom algorithm is much faster than the 
traditional relaxation-based technique. The number of 


propagations in every polygon p; 1s reduced from A; (the area 


of p;) to A;/(w,**2), where w; 1s the weight of a unit distance 
71 els Cheah 

The user can change the threshold parameter in straight- 
forward algorithm, as he can also alter the stride parameter, 
in addition to the threshold parameter, in the modified 
algorithm. The threshold parameter ensures that the algorithm 
stops as it 1S with some threshold for the optimum. The stride 
parameter dictates the minimum weight of an edge thus making 
bigger strides in low cost areas. The stride parameter 
indirectly controls the number of edges and nodes to be 
eliminated, hence making the algorithm faster. 

An adventageous side-effect 1s the elimination of some 
nodes. These nodes are eliminated if they are not connected 
to more than one node (both source and goal points cannot be 
eliminated). The psuedocode for these procedures are given 
in the same order as in the program: Horizantal search, edge 
eliminating, vertical search, edge eliminating and node 
eliminating in Figure 14. Furthermore, a continuation to the 


example in Figure 12 is given in Figure 15. 
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Algorithm: To find the Weighted Region Least Cost Path 
in the modified version of relaxation-based 
Grid alqorteam 


Given: Two dimensional grid of points: GRID 
SOUrTGEeE DOInt seo eiS) 
Goal point: (1g,j9g) 
TAGesn@ ld ak 
Se eiela* 


Queput: Cost, of minimization Erom Goure=e to coal) eos 
The minimum path Erom Sotimee tougoal. 
(1Syg S) (ver ee: 


procedure ; 


--This code 1s inserted at point xx in Figure 11 

Horizantal Search and Elimination(HSE) 

Vertical Search and Elimination(similar to HSE) 

Eliminate the nodes which are not connected to more 
than one neighbor 

--The rest of the algorithm is the same except 

--while loop uses procedure propagation to propagate 





Figure 14: The Pseudocode for Modified Grid Algorithm 


C. RESULTS 

We performed a series of experiments using the Meridian 
Ada compiler in Sparc stations in order to realize how gene 
parameters effect time and cost. Some of results from these 
experiments are combined and shown in the Table I-IV. In Table 
I and II we tested threshold parameter and concluded that the 
greater the threshold value is, the faster the program is. 
Changing the threshold value does not have a significant 
effect on the least cost. The results taken from the 
experiments with modified program are presented in Table III- 


IV. In these experiments, we kept threshold value constant 
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(equal to zero) and changed the stride value between 1 and 
moo. Im Gee first row we realized. that thes cost of 
modifications increased the amount of time for computations. 
As seen form the last column, the time gets less even though 
the second trial took longer time. The number of the nodes 
increases depending on the stride value. Meanwhile, the 
disadvantage is that the least cost gets higher for higher 
stride values. 
D. DRAWBACKS 

The modified algorithm requires more overhead. The over- 
head as reflected by the results 1S worthwhile because the 
overall running time of the modified algorithm is smaller 
Eee the Graginal algorithm for an approximately close 
stride value. Time overhead might be reduced by not requir- 
ing that distances between all neighboring gridpoint be pre- 
@Senoueea,. Again, this algorithm is intrinsically suboptimal. 
There 1s a possibility that the graph is decomposed into more 
than one connected component due to the elimination of edges 
and nodes. If such is the case, then source and goal points 
might not be reachable from one another. This problem will 
manifest itself by having the minimum cost of the grid not 
changing from the assigned initial value. If there 1S more 
than one connected component, and the source and goal points 
are in the same connected component, then there should be a 
technique that avoids computations in the other connected 
components. 

The program codes for the modified grid algorithm are 


enclosed in Appendix B. 


Zo 


STEP 1. Horizontal search is performed 


(e.g. with stride = 10, in this case 


the weights an edges are added until 


1t gets more than stride.) 


STEP 2. Edge Elimination 1s performed 





Figure 15.a: An Example For A Modified Algorithm (T) 


IAS 


ROURCE 


(COnr nuee a taom br gure 35.4) 
STEP 3. Vertical Search is performed the same as 


MOm@uZoncal, search 


STEP 4. Edge eliminating is performed. 





Figure 15.b: An Example For A Modified Algorithm (IT) 


LG) 





(Continued: frome Figquvem rae: 
STEP 5. Node Elimination 


STEP 6. Cost Minimization and ftingune 


least cost path. 


Figure 15.c: An Example For A Modified Algorithm (IIT) 
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Table 1: EXPERIMENT SHOWING THE COST AND COMPUTATION TIME 
VERSUS THE THRESHOLD WITH CLASSICAL ALGORITHM FROM (1,1) TO 
(80,80) 


ee aoin— 














Table 2: EXPERIMENT SHOWING THE COST AND COMPUTATION TIME 
VERSUS THE THRESHOLD WITH CLASSICAL ALGORITHM FROM (10,1) 
TO (55,80) 


ee a 
CO 













DS, 


Table 3: EXPERIMENT SHOWING THE COST AND COMPUTATION TIME 
VERSUS THE STRIDE WITH MODIFIED ALGORITHM FROM (1,1) TO (80,80) 


Number of 
Threshold Stride Cost Time(sec) Nodes 
Deleted 


CN 
CC Es 
CE 
CE 
of 6 as ore 


Table 4: EXPERIMENT SHOWING THE COST AND COMPUTATION TIME 
VERSUS THE STRIDE WITH MODIFIED ALGORITHM FROM (10,1) TO 


(55,80) 
Nodes 


Threshold Stride Cost Time(sec) 
Deleted 


CC OC 
CS CS EC a 
CC EE 
CS EEC 
CE 
a CE 


* = QOvernead Causes ENiSs ine reace ene mn. 










Number of 
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IV. THE ENVIRONMENT 


In this chapter we describe both the hardware and the soft- 
ware environments in which we implemented the parallel algo- 
rithm presented in Chapter V. The first section introduces the 
Transputer [TRA87], and the second section introduces Alsys 
Ada [ALS90]. A guide for program development using Ada on the 
Transputer 1S presented in the third section. 

A. TRANSPUTER 

The Transputer implementation 1s based on the concept of 
the Communicating Sequential Processes(CSP). In order to uti- 
lize the transputer effectively, we need to understand the way 
it works. Transputer 1s a microprocessor with its own local 
memory storage and four links designed to communicate directly 
with other Transputers. The larger the number of processors 
in the network 1s, the more processing power, the more memory 
and links are available. The difficulties also grow with the 
number of processors. The most visible difficulty in the net- 
work 1s to avoid deadlock in which communication fails and 
results in processes walting forever. 

There are different types of Transputer: T2 (T212, T222), 
hoe ee 1425) oN Te (every t301, T8005). We worked with T800 
Transputers. A block diagram of a T800 Transputer 1s presented 
7 Figure 15. The major components of T800 Transputer, as 
seen, are memory, processor and communication system connect- 
SGue tamer 2 be Husk 


The high level programming language OCCAM [0CC88] [0CC89] 


Sh 


is the primary language used for programming the Transputers. 
It is designed to run concurrent processes on a network of 
Transputers. Concurrence and communication are two main con- 
cepts in OCCAM. They allow processes to run simultaneously and 
transfer information through channels from process to pro- 
cess. Processes communicate by message passing, do not share 
variables, and synchronize only when they communicate. Commu- 
nication 1s synchronous and unbuffered. 

The host computer for the Transputer network that we use 1s a 


PG 208 


FLOATING POINT UNIT 


CED PROCESSOR 


I 

N 
T 
E 
R 
N 
A 
1 
B 
U 
S 


alae EXTERNAL MEMORY 
INTERFACE 


Figure 16: A Block Diagram of T800Q. 
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B. ALSYS ADA 

Alsys Ada Compilation System, which consists of a compiler 
ae Seo iedeq mms sUSeQutem handling Ada programs in a 
Transputer programming environment. We used Alsys Ada as the 
language of choice for developing our WRLCP problem code. 
Below we describe the most salient feature that distinguishes 
Alsys Ada usage from ordinary programming in Ada. 

1. Channels 

Communication between Ada programs 1s provided by using 

transputer channels via the implementation defined package 
CHANNELS. The CHANNELS package contains a generic package 
CHANNEL_IO which defines input-output for values of a 
specified object type. READ, WRITE, READ OR_FAIL, 
WRITE_OR_FAIL procedures within this generic package are used 
for input and output between channels. The distributed 
application 1S written as a set of independent programs for 
Single or multiple Transputers and communicate through 
channels with unique names. 

We write a COMMON package, which contains declarations 
common to more than one Ada program and which also contains 
an instantiating of CHANNEL_IO to allow data to be 
communicated between independent programs. The data type of 
the channel, common to an application, is defined in a COMMON 
package. For example, to declare a channel of new NAMEI1 of 
type DATA TYPE, we include the following in the COMMON 


package: 


a 


DATA_TYPE : (can be any type: integer, records.) - 
package NAME] is new CHANNELS.CHANNEL_IO (DATA_TYPE) ; 
In the programs it is used as follows: 

declaration-- Ada program segment begins here with declaration 
NAMEZ2 : CHANNELS .CHANNEL_REF : =CHANNELS .IN_PARAMETERS (vi 


rtual channel number); 

--this channel is used for input trom the other processors and virtual 
channel ---number ie given ly the programmer(eé.g. i or 2 ae 
NAME3 : CHANNELS .CHANNEL_REF : =CHANNELS.OUT_PARAMETERS (virtual 
channel number); -=this is for output to the otter se eeessens 

begin ~-main program starts here 

NAME1.READ(NAME2, the same data_type) ;--get the input from the 
thannel --NAME2 

NAME] .WRITE(NAME3,the same data_type);-- put the data_type 
value to --the channel NAME3. 

end 

2. Harnesses 

In order to run Ada programs in parallel on a single 
processor or a multi-transputer network, we need to use an 
interface, which 1s an occam process called a harness. A 
Nnharness 1s used as a wrapping for the Ada program to be 
accepted by a Transputer. For every Ada program, two occam 
harnesses have to be created. Harnesses are explained in more 
details in section C. 

C. GUIDE FOR PROGRAM DEVELOPMENT 
In this section we present some helpful points to make 
program development easier. After learning the MAKE Program 


Maintenance Utility, the tools for checking the network, and 
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studying the examples, one will be ready for program 
development in this environment. 
1. MAKE program maintenance Utility 
MAKE 1s a utility program designed to help control of 
programming environment, EO automate the process by 
determining which parts of the program is changed since the 
Hasctmecemellatzion and rebuilds them accordingly. 
makefile 1S a script file, written by programmer and 
directs MAKE. You can find MAKE commands below: 
make family: Creates the Ada family and library sub- 
directories. This 1S a one-time-only operation. 
make: The standard command for building the executable 
programs after changes have been made to the source. 
make run: Executes the compiled program. 
make help: Displays the MAKE commands. 
make -n: Displays but do not execute commands. 
make check: Checks transputer topology. 
make clean: Deletes all files except source files. 
make *.0: Make Ada object codes. 
There 1s a batch file named doit_all that executes the 
first three commands:make family,make, and make run. 
2. Checking Tools 
Besides make check there are some more executable files 
which check the network topology: worm.exe, chknet.exe. 
3. Examples in Steps of Instructions 
Tt 1s better to start with complete examples to get used 


Eo te 


ct. 


I. Make your own directory “and eopy @sne — een_rice 
installatten mee Ve. 

>copy d:\alsys037\source\generic\*.* 

Complete documentation can be found in the files 
read.me and show.me. 

II. In this environment, there 1s an Ada source file 
proj.ada containing procedure proj. If you edit your own code 
with these names, 1t means you are ready to compile your code 
in Alsys Ada environment. 

III. You can type doit_all which executes three 
commands: 

make family,make,make run. 

Iv. It ais time to try the examples on a single 
transputer and then on multiple transputers. You can refer to 
the appendix C and Alsys Ada User Manuals. 

V. Communicating Ada processes on a single transputer 
needs the list of files below: 

makefile: A script file written by the programmer and 
executes the commands according to makefile. 

family.inv: This file creates the library environment 
(It does not change.). 

proj.inv: This file directs compiling andeemmaing. 

main.oce: In order to integrate Ada with other 
languages a well defined interface is required. Ada programs 
may then be run in parallel on a single processor or 
distributed across a multi-transputer network, just as occam 
processes. This is a default occam harness provided as part 


of the compilation system in both source and compiled forms. 
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The main body of the harness consists of three processes 
operating in parallel: 

- A multiplexor which combines the error output and the 
Standard output of the Ada program. 

—— (eoeeemmehannel Collector which colleces any output 
from the error stream and routes it to the standard output 
stream of the server via the multiplexor. 

- A process which sets up the input and the output 
channel vectors of the Ada program and then invokes it, 
informing the other processes upon completion. 

merger.oce: This default harness is used to collect the 
error output from up to some number of Ada programs and send 
1t to the standard output stream of the server (It does not 
change.). 

projh.oce: Each Ada (here PROJ.ADA) program has 1ts own 
mini harness which provides a clean interface to the program 
in terms of the channels used. Main harness 1s used to invoke 
each of the mini harnesses in parallel. 

projh2.oce: This 1s the dummy harness required to allow 
linking of a foreign Ada program with the occam libraries. 

main.ink: Gives the file list to link. 

VI. The files needed for multiple transputers are 
mainly the same but they should be modified according to the 
network and presented below: 

makefile 

family.inv 

proj?.inv: They should be as many as the number of 


transputers. 
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mainh.occ 

merger.occ 

proj?h.oce: They should be as many as the number of 
Cransputeers., 

proj?h2.oce: They should be as many as the number of 
transputers. 

main.pgm: This is the only file not needed for one 
Single transputer. It describes which virtual channels are 
equal to which physical channels. 

main.lnk 

An example is provided in Appendix C containing all 
these files, and it does not need to be changed for similar 


applications. In Figure 16 all these relations are shown. 


aes: 


© 
=: 
= 
O 


test_fam ; 
we bind 


bin 


common.ada 


ompile 


rojl ada proj0.o 


nk 
proj lh.c&s 
proj lh.m8s 


compile 


rojO.ada 


a 


ili 


proy lh.oce occam 
proj lh.t&s 


projlh2.oce [Occ proj 1h2.tax 


jlink icon 
occam iserve 


projOh.occ occam 


projOh2.occe PCC projOh2.tax 


occam 


mainh.oce mainh.t&s 


Figure 17: Diagram Of The Steps Involved In Program Development 
Using Alsys Ada On Transputers 





be 


V. A MULTICOMPUTER ALGORITHM 


This chapter introduces a parallel algorithm for the WRLCP 
problem. In the first section we present the implementation 
on an INMOS Transputer network using Alsys Ada. In the second 
section we discuss a variant of that algorithm that uses less 
Communit Cat On Veter ser 
A. IMPLEMENTATION 

To simplify development of the algorithm, we made every 
possible effort to have it treat all processors symmetrically. 
This approach allows the algorithm to be scalable for a 
different number of processors without much change. We show 
the pseudocode for a version of the algorithm running on two 
PYrOCeESSOrs 1m JFPigure. 7. 

The Root Transputer, which is the only Transputer having 
direct connection with the host PC, reads the data file and 
sends equal portions of the grid to other processors. At the 
beginning of the computation the Root Transputer sends the 
values of the threshold, the stride, the number of processors, 
the source and goal points to every processor. 4) eae 
processors generate the weights on edges. All of them start 
cost minimization at the same time. Every processor updates 
the values of the grid points in its portion Of the orld ee 
the points that are at the border of the grid portion assigned 
to a processor, the costs of the data obtained from the 


neighboring processor at the previous iteration are used. At 
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Blcqonrtenm.. ©o find the Weighted Region Least Cost Path on 
Cwo processors 
Given: Two dimensional grid of points: GRID 
DOlmeer pointe as, 7s) 
TOSI clon ac ame mare h 
Wilcesiold. . 


Mucput. COst Of minimization from source to goal: Cost 
The minimum path from source to goal: 
(is,js) ine 1) 


procedure: 


Initialize 
Read the data file 
write(initial grid data) 


for all nodes loop 
SGalculate weight (1,])} 
end loop 
write (source, goal,n,p, threshold) 
Prag = true 
while flag loop 
flag = false 
for all nodes loop 
write (border costs) --exchanging the costs on border 
read (border costs) --nodes 
find minimum value 
change = node(1,j).old_cost-node(i,j).current_cost 
if change > T then 
raise flag 
end if 


end loop 
write(flag) --checking termination code 
read (flag) 
end loop 
read(minimized costs) --for integration the solution, it 
--gets what the other processor did 
output the WRLCP 





Figure 18: Psuedocode For ebact Je plore a 


From The One Processor”s Point Of Vi 
The Other Processor aa Have Read ae Writes 


Interchanged. 


the end of the computation of an iteration, neighboring 
processors exchange values of the border grid points to be 


used in a later iteration. After every iteration, the 
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processors, as a whole check each other to decide whether to 
stop or to continue. If they decide that an adequate solution 
has been found they stop. All the processors write the 
obtained costs back to the Root Transputer which then displays 
the least cost path from the source €o the goal points Jenene 

The pattern for “WRITE and READ” between the processors is 
shown in Figure 18. Other correct patterns exist, but it 
should be emphasized that a very important issue is deadlock 


avoidance. 


PROCESSOR-2 


INITIAL GRID DATA 


—— 
SOURCE, GOAL, N72 THRE SsHonm 
WRITE 
a 
EXCHANGING BORDER COSTS 
READ WRITE 


ee aes 
FLAG 


Sencar aon READ 
HECKING TERMINATION CODE 
READ —.. 


SS MINIMUM COST VALUES LS 


Boas at 
INTEGRATING SOLUTION 


Figure 19: Write and Read Patterns 
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We have a running program on two processors, which is 
presented in appendix C. Furthermore, the methodology can be 
applied to more processors. Different patterns can be used to 
Basta elOnmmeme Gridegrapm on a numberwof processors.in Figure 
19 it is shown how the grid can be partitioned in different 


possible patterns, e.g. for four processors. 


processors processors processors 





Figure 20: Patterns for Partitioning 


B. A VARIANT ALGORITHM THAT SAVES ON COMMUNICATION 

The cost of communication can be very large since we have 
to exchange all the border values and checking parameters 
between the neighbor processors throughout every iteration of 
cost minimization. In reality, a lot of values exchanged 
across the border of the grid portions are redundant because 
they can remain unchanged through more than one iteration. 
That is why saving on communication becomes very important. 
As we try to speed up the algorithm, the cost of communication 
SNeGlMlGmemor "Slow te Gown, especially if part of the 


communication is pragmatically useless. We present a modified 
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algorithm that takes advantage of the previous observation. 
In the modified approach, these border points which we changed 
in an iteration are marked. At the beginning of the data 
exchange between neighboring processors, each processor 
informs its neighbor about the number of the border points 
that changed. Then it proceeds to send the low cost values for 
only these points. More saving on communication can be 
achieved by using variant records. Unfortunately, we could not 
set this to work in the current development environment .The 
pseudocode for this approach is presented in Figure 20. 


Algorithm: The variant algorithm that saves on 
communication 


procedure: 


for all border points node(i,j) loop 
change = node(i,j) .cld_ cost = node(i4) currence ucece 
if change /= 0 then 
mark node(i,j) 
end if 
end loop 


count = number of marked nodes(1,j) 
write (count) 
for 1..count Toop 
write marked node(i,j) 
end loop 
read (count) 


for 1..count Woop 
read marked node(i,j) 
end loop 





Figure 21: The Pseudocode For The Communication-Saver 
Variant Algorithm Of The Parallel Algorithm 
From The One Processor’s Point Of View. 
The Other Processor Will Have Read and Writes 
Interchanged. 
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VI.CONCLUSIONS AND RECOMENDATIONS FOR FUTURE 
RESEARCH 


We have developed an efficient version of a grid based 
algorithm to solve the WRLCP problem. Our algorithm shows a 
Significant decrease in the computation time in comparison to 
the original algorithm. The experienced loss in solution 
pi@@lieacy ptommoOl Proportional tO the saving in computation 
time. 

As a step towards developing a parallel version of this 
algorithm on a network of Transputers using Alsys Ada, we 
started developing simple parallel grid algorithm for the 
WRLCP problem in the above environment. Due to difficulties 
in dealing with that environment and due to lack of time, we 
stopped at the stage of developing adequate parallel 
algorithms for that environment, hoping that others will 
pursue our efforts towards the initial goal. 

The main emphasis in our parallel algorithm was compile 
time partitioning and mapping of data. Garcia ([GAR89] 
implemented the parallel algorithm (Local, Asycnhronous and 
Iterative Parallel Procedures (LAIPP) Algorithm) presented in 
[SMI88] on a network of Transputers using Logical C. 
Scheduling was dynamic by farming out computations to 
avallable processors. Garcia’s results showed that at a 
certain point, increasing the number of processors decreased 


the speedup. We attribute this to excessive communication 


A5 


delays involved in the scheduling. We directed our effects 
towards nearest-neighbor patterns of communication, and we 
believe that this appropriate approach to handle this problem. 

Many problems need yet to be solved. As a starting point, 
software tools for automatic generation of Ada harnesses, and 
automatic mapping of Ada programs need to be acquired. This, 
and upgrade in the existing hardware setup will bring about 
more rapid program development. Currently, the process of 
program developing in that environment 1s extremely tedious. 

First order improvement to the existing parallel 
algorithms can be attained by usSing variant records for 
communication across channels. 

More serious improvement include: 

- designing an asynchronous version of the parallel 
algorithm (data will be communicated only when needed), 

- uSing queues or heaps can decrease the computations, 

- uSing the routing library developed in [FAL92] might 
provide us a way to compare our algorithm to more efficient 
algorithms that are not constrained to nearest neighbors. 

During the course of our work, we encountered problem with 
theoretical flavor which yet to be solved: 

- ensure that edge and node elimination in the algorithm 
(to reach a parallel analog of the modified sequential 
algorithm) does not separate the grid graph into different 


connected components. 
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APPEND 


--Title 

== le lavone 
--Date 
--Revised 
==(Ouls ec 
==COMmMpn Herm 


IX A: SEQUENTIAL ALGORITHM SOURCE CODE 


STRAIGHT-FORWARD ALGORITHM 
CENGIZ EKIN 

2UeaviG7 92 

04/10/92 

THSSits 

MERIDIAN ADA 


--Description:Reads data from file,input starting and goal 


per Ares 


WLlth  “TEXT28O, OS. 1 ven e- 
use TEXT. boo Saher 


procedure MAINI 1s 


--finds the minimum cost path. 
TASK “CONTROL, CALENDAR 
CALENDAR ; 


package INTEGER_INOUT is new INTEGER_IO(INTEGER) ; 
package FLOAT_INOUT is new FLOAT_IO(FLOAT) ; 
use FLOAT_INOUT, INTEGER_INOUT; 


START TiMiEa 
END_TIME 


Sp meh Gn CP. ra Galle | 

Q 

VOLTA, VOLT 

NOP 

EA COUNTEE 

EA oo 

type ELER is 
record 
WEIGHT 
DISTANCE 
end record; 


type GRID_POINT 


record 


CURRENT eCOs il. 


OLDECest 
N, 
E, 


INTEGER; 
INTEGER; 


FLOAT ; 


> INTEGER: 
SR ENG (ee a eee 
>: INTEGER; 

7 INTEGER = —oU- 
[INTEGERS 


= SUPE CI, 3) 


INTEGER; 


30 


Sy 
W qe ee 
Pao bon 2) SHE SGU: G GliRenp sy ies 
end record; 
Spe vce Pers array (i. (NOP +1),0..(NOP+1)) of GRID POINT; 
B Chip: 
NE eee Ober Ee. 


moe oPeEAGr ts array (1..1000) of INTEGER; 
Byoe oer Gril is array (1..-1000) of STRING(1..5); 
type STACK 1s 


record 
Soler So LOR AGE > 
LATEST Vr ECER: 2= 0. 

end record; 

type STACKI1 is 

mee eed 
STORE JoLORAGE 
ee Eh J PNTECHER es = 10 


end record; 


S Sars JENG Grol eS J NG ailie 


eo 


CO 


Mmagcedtrewruon (S >: 1m out STACK; EF : in INTEGER) is 


begin 
Serna f = o>. LATRol + 1; 
SPolOREVogbAT Roly =:= E> 
end PUSH; 


mmm www ww i ie 


muee@eauie PUSH! (Sl- in out STACK]; El : in STRING) is 


begin 
Sie eAthoee c= cGleLATEST + L: 
Sl STORE (sl LATEST) <= El; 
end PUSH1; 


me ae ee ee SS Pe ee 


— mmm mim ss S| 


oul 


procedure POP (S : 1n out STACK, 9H Somes ethene. 


begin 

E += S.STORE(S2EATEST)= 
S.LATEST ;= S LASES Ly 
end POP; 


procedure POP] (Sl 3 in cut STACK]: Ei veers: fie 


begin 

alee sh SOROS CS Less) 
Sl. LATEST: 2]="Si2 tates ia 
end POP1; 


-- This function computes the execution time (CPU TIME). 
Function CLOK return BEOnt as. 
Funct On VCLOGCK. Bet ir me semen 
pragma INTERFACE(C, CLOCK); 
Is oo one 
> : LOA 
begin 
task_control.pre_emption_of €; 
T; = "CReecr; 
1 f ES ene 
raise TIME_ERROR; 
else 
S :=-FLOAT(D)/ 1. Cac. 
end if; 
task_control.pre_emption_on; 
returns 
ena we h@nK- 


meee eee sss es es eS 


procedure CAL_WEIGHT (I,J : in INTEGER) is 


begin 
tf i] see then 
Bil,d) N2WweEtGe tai. 
Bil -d.) NSP LS tANGE fake 
B(l+i,7). CURRENDSCOST. = 10000. 
else 


B(I,J) .N.WEIGHT ;:=1+ ABS(B(2, 40) 3eUnt 2622. 


5) 


Siem eURRENT COST) ; 


end if; 
twee = NOP then 
Pel) 5. SeuGHE = —1; 
Bie) DiS teiNGE..= —l]> 
Peele de CURRENT COST e=:=10000 ; 
else 
Peo. Vetehe  s=) leABS(B(1,J}) .<CURRENT_ COST - 
Pet wt ee CURRENT COs.) 
end if; 
Pea listener 
Baked) ss. WELGHT <= => 
Batted Sb ko LANGE = =]: 
bit Oe oUKne it scOosSr 2=LO0000 > 
else 


Pew vom eens sl sab Steir, J) -<CURRENT_COST - BiI- 
mr SCuRRE NT 2ceST ) ; 


end 1f; 
fio = |. then 
Pew WeN EGET 2S.-=1- 
Pier) ev DESTANCE = => 
Biverd —WimeukneNnTr COST <—0000- 
else 


Peer Ve rGh hee =i ABS eww) “CURRENT COST — B(L,J- 
ie CURRENT COST) ; 
end if; 
end CAL WEIGHT; 
pasececume FIND MIN (I,J : in INTEGER) is 
begin 
Peme bt ) CURRENT COST*s> fe (eed) SCURRENT COST 
+abs(B(I,J) .N.WEIGHT)) then 
(Srl CURRENT SGO@sT 
+abs(B(1I,J) .N.WEIGHT) ) ; 
ON GSI sie 


eo 0 ORM Olean COs ag 


Bio veo hence lr TON 
SyqG, aloes 
fee Bil) CURRENT COST > C2 ie CURRE NL eGOs tT 
+abs(B(I,J).E.WEIGHT) ) then 
BIE tL eS Ofededay nt Ce Suk 
tabs (B(1I,J) .E.WEIGHT) ); 
qlab sed Gece 


Pie eeurRENir COST 


El ee > Reel ON 
end if; 
Zeiss CURRENT COST =& eee — ee) CURRE ND 2eoOsT 


I 
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+abs (B(I,J).S.WEIGHT)) then 
(B(E-ly) CURrE Rae Osi 

+absi BL Mees Velen eb) 

Y SCURn 


Bil, d) -CURRE ees 


Bild) 2. DiEREC hae 
end if; 
1£ B(1L,J) -CURRENZSGes ia (Bae —)) SeuRR Ee aeocar 
+abs(B(1I,J).W.WEIGHT)) then 


I 


B(I,J) .CURRENT_COST := “(Bigi( 3 CURR Ee eum 
+abs(B(I,J) .W.WEILGHT) ); 
B( Ld) DEREC REG 7a WES ae 
end 1 


end FIND_MIN; 


TT 


cece eee i 


while © = 47" sleep 


COUNTER: ake 
for fan ee Neer lcer 
fOr clean NOP loop 


OLS LANCE = .- 


i 
( N 
2E. Dist Aner. 
S 
W 


Beh.) 
Bee) 
Bel 
Bie) 
end loop; 

end loop; 

OPEN (INF ,MODE =SeIN FILE, NAME =— “temas 
PUT ("THE DIMENSION OF MATRIX ne 


(DES TANG ha =4 
is tAN Caw — a 


GEL, ere 

PUT_LINE ("ENTER THE VOLTA (the optimizataen tolerances 
PUT (“VOLTA =") GE Olive 

PUT_LINE ("ENTER THE SOURCE POINT !"); 

PUT ("Sx = ) >GET 60S) PUT (2 SY U= 2 GE lca 

PUT_LINE ("ENTER “THE GOAL POLNT 773 

PUT( "GX = ") 7GBR (GX PUT ( (CY o-eee Chl (eae 


—_——_—— eee ss ea 


-- CLOK function begins to compute the execution time. 
STARTS Eine = Chere: 


— eee wwe eee es ee es 


—_ tee eee eee ee eee eee ss Ps es lr: ee 


for ROW in 1..NOP loop 
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Eomecemheran ise.NOP loop 
Chee clNe. B( ROW, COL) .CURRENT COST) ; 
end loop; 
end loop; 
--CLOK function finishes the computation of execution time. 
Pia eEME 2:= CLOK: 
Piero Atb TIME TO READ THE DATA FILE IS : ue. 
PUT (END_TIME-START_TIME,4,4,0) ;new_line; 
Sore EMP <= CLOK: 
--It determines the borders and calculates the weights of the 
edges... 
ferret ean 1..NOP loop 
Bowesahn  l..NOP loop 
CAL WEEGHT(1L,J) ; 
end loop; 
ee, Loop; 
--It makes the costs max number in order to use them in 
comparisons for finding 
-- the minimum.... 


fom rt un 1 NOP loop 
Bogue ines. NOP’ tlLoop 
Piero CUR nem scOST. == 10000); 
Ve Onbmee@ot = sas il wd) CURRENT COST: 
end loop; 
end loop; 


B(SX,SY) .CURRENT_COST := 0; 
while COUNTER > 0 loop 
COUNTER := 0; 


Pomel eat La.NOP loop 
Hor weih eNOr Loop 
PND mehr aa 
Otel esol COSne eid) CURRENT COST; 
tim Owes VOUrA then 
COUNMERe. = COUNTER +h. 
end 1 f; 


DD 


B(I,J) .OLD_COST := Billy) 2CuUrRE ize e..u 
end Loep, 
end. keep; 
end Loot: 


Se eee eee 


leon 
PUSH (S,B(GX%, GY) GURREIT Goer] 
Pura it 3 GY) .-S LR eCLION) = 


( 
SQ, =B4G Y) sPreBerrel 


1£ SO oe then 
PUSH(S; B(GX, GY} Na TSTANGE 
C6 — eC eo 
elsif SO = “BAST sinew 
PUSH(S,B(GA%, GY) 25. Distaner 
GY “t= 4G yee, 
elsif -SO’ = "SOUTH" “Ehen 
PUSH(S; BUGA, GY )o%52 DD lbSolANee 
GX := GX-1; 
elsit SO = “WEST wemen 
PUSH (S, B{GX, GY) (Webi StTANGeEe); 
GY 32] 4Gye. 
else 
exe. 
Syavel, see 
exit when GX =SxXsandsGy = oe 
end loop; 
PUT SEN Ey DISTANCE CoOsrT DIRECTION ==: 
PU a Ee nn Sa ss 
loop 
POPS EB) e2UT vb ese hee a ae 
POP (SE) PUT (Ee). Un we Be) 
POP CS abn: 
Di EL = “NORTH weno 
PUT ("SOUTH") ;new_line; 
elsif El = "EAST “ then 
PUT ("WEST") ;new_line; 
elsif ElG=3-7 SOUDRH" Semen 
PUT ("NORTH") ;new_line; 
elsif B= “WEST “semen 


PUT ("EAST") ;new_line; 
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yinvel abie % 


Sexat whem S.LATEST =0; 
end loop; 
ShOSE (INE) - 
PSV weiter Pon Limisnes the Computation of execution time. 
Bie Peo = CLOK: 
Poe Tora PRIME TO EXECUTE THE PROGRAM IS : gd 
PUT (END_TIME-START_TIME,4,4,0) ;new_line; 


Powel eoOnk MORE TEME -:y(eS) or no) :");GET(Q) ; 
Pa Oe- on “phen exit: 

end if; 

end loop; 

end MAINI; 


arg 


APPENDIX B: MODIFIED ALGORITHM SOURCE CODE 


--Title : MODIFIED ALGORITHM 
--Author : CENGIZ EKIN 

--Date 1020/7 OO7 OZ 

--Revised : 04/10/92 

== eUurSse § 2 Siar sigs 


--Compiler: MERIDIAN ADA 

--Description:modified program, searches and deletes nodes 
and ----edges 

with  TEXT_1I0,0S_ TYPES, TASK. CONTROL ee ore lieAk 

use TEXT <10,0S._ TYPES, CALBENDARS 


procedure MAIN2 1s 
package INTEGER_INOUT 1s new INTEGER_IO(INTEGER) ; 
package FLOAT_INOUT is new FLOAT_IO(FLOAT) ; 
use FLOAT_INOUT, INTEGER_INOUT; 
SAR at i Mirae 


BND? EME * sh LOA. 


SX, SY,G4,GY, 2)/CO2, COs, C04, COs 2 i TEGEr 


BE, COUNTER nw. 5 » INGEGER. 
Ede STRING! oe 
QO ‘STRING UI ie =e 
DERI VOrmAy Olek > INTEGER; 
NOP INTEGER =. 007-5 0 
MARK : BOOLEAN; 
type ELER is 
record 


WEIGHT : INTEGER; 
DISTANCE = INTEGER: 
end record; 
EvVpe GRIDSPOLNT was 
record 
CURRENT 2C@si- 
OLDseOs7 : INTEGER: 


Bie 


ba] 


=) el = 


Pie: 
Pleeer ron meotRUNG(] 25) 
ACTIVE a) a COLORS AN 


end record; 
oer lemtema- rea 0... (NOP +1)70..(NOP41)) of GRID POINT; 
B oe Le 
TNE 2 Loe TYPE 


type STORAGE is array (1..1000) of INTEGER; 
type STORAGE] is array (1..1000) of STRING(1..5); 
type STACK 1s 


Becord 
SORE oO LORAGE: 
BAS P NE Gr Ree 0): 

end record; 

type STACKI1 1s 

record 
SATO} >: STORAGE] ; 
LATEST TUNE GE Re: = Or 


end record; 


S eolaek > | : STACK: 


www www se 


procedure PUSH (S : in out STACK; E : in INTEGER) is 


begin 
Se oewesoeer—= SG. LATEST + 1; 
SV olORe (oebAlTEGLr).:= E> 
end PUSH; 


—_— eee eee ees SEs S| 


Breocedure PUSH! (Sl : in out STACK]; El : in STRING) is 


begin 
Sieeoieot. -= SIPELATEST + 1; 
Sl. SEROUSGE (St bye Mau ae oe ee 
ena PUSH: 


2) 


wm meeeweei ii i  ee i le Ce 


procedure POP (5 : 1n OUE STACK 3s). outs ©) Reem ny me 
begin 


E <= S.STORE(S.LATEST)-: 
5. LATEST <= SLATES 
end POP; 


procedure POP] (SH: in Out STACK] = Bec ieee ee 


begin 

Bl := S1.STORE( > Ulla. 
Sil; LATEST 3] Sl yEbAtbSr aaa 
end POP1; 


eee i i eee 


-- This function computes the execution time (CPU TIME). 


function CLOK TreturneeLesT srs 
LUNCEIOn ChOCK serum. ine 
pragma INTERFACE(C, CLOCK); 


IBIS Fs oker 
S - FLOAT 
begin 


task_control.pre_emption_of f; 
=. COC ke 


LE T = =) then 
raise TIME_ERROR; 


else 
Seek (Tee ee Eo. 
end if; 


task Gene eol= present loneaon: 
EFSlULaes. 
end LOR; 


cm mcm mmm mmm eee eee ees ES es 


procedure CAL_WEIGHT (I,J : in INTEGER) is 
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begin 
if ieee) NOP then 


boo Ne weGaT == —] > 

Eee hie bomen bh = = 1 ; 

ote eeOleane iro COst :=10000;> 
else 


See) a ene: — lee boi Bir, ) CURRENT COST - 
Eaal+ il, 30) .CURRENT COST) ; 


ends i: 
ite = NOP ehen 
Bethe ee awe lGh fe. — 1) 
Ea eee eon ANCE os =]: 
Sole) seu eeNr COST 2-=10000- 
else 
Stew See GH hes P+ ABS (B{I,J).CURRENT_COST - 
Poe !) CURRENT COST); 
end 1f; 
Ie po I elateyal 
Eel ormomVe ren hs t= -1- 
El eo eS RaNGe, <=) a): 
Bul -iaw) CURRENImeeSsr --10000- 
else 


B(I,J).S.WEIGHT. :=1+ ABS(B(I,J) .CURRENT_COST - B(I- 
iw eee URRENT COST) ; 


end if; 
ihe then 

Bw Wo VEE Gea = 1 

Piel a oie LaNC Bae — 1: 

Peel) - ee URBENT COST :=10000; 
else 


Peo NE eee ie. = ees (1 ee eRRENT COST - Bi(1,J- 
Po OURRE Nr COST); 

eng. if; 
end CAL WEIGHT; 


ewww wwe wee 


procedure FIND_MIN (I,J : in INTEGER) 1s 


begin 
1f B(I,J) .ACTIVE = TRUE then 
Pees eed Ne pietance /= -1 then 
Veet J) VeURRENT COST > 
Tee fees el a) ie DES PANCE) ) J) .CURRENT_COST 


+abs(B(I,J).N.WEIGHT)) then 
Eee CURR eNT COST == 
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(B( (I+abs (B(1,J) .N.DISTANCE))) 73). GURREI|iaee un 
+abs(B(I,J) .N.WEIGHT) ) ; 
B(i,d) DIREGIToe)) = S NOR wee = 
end if; 
end if; 
Lf Bli,Jd) . Boise — ene 
Lf Blt, 3) -CUuURRE Ie Ceo i 
(B(L, (J+abs(B(1,J) EB. BErSTaANCE) )) -CURRE taee om 
+abs(B(1I,J) .E.WEIGHT)) thea 
Ba Lj) ) PCuURRE NDT AeG@ oc .— 


(B(1L, (J+abs (B(1I,J) .E. DISTANCE) [euUR EE eo a 
+abos (3B (iad) eeWwEiGHh Th): 
B( il ,J) .DIEREGZ LON = a ee es 
end if; 
end if; 


1£ B(I,J).S.DISTANCE /= -1 then 
LEB (isd) SCURREN ies ihe (B((1I- 
abs (B(21,J) -<S. DISTANCE) 33) ee eel ace sy 
+abs(B(I,J).S.WEIGHT)) then 
B(L,J). CURRENT COST... 2iee0— 
abs (B(1,J).S .DISTANCE)) 7U) SCURREI i =eo = 7) 
tabs (B(i,)) -S2yEten ris 
Bil, J} DIRECTO =. oO Ula 
end if; 
Sygiel 3. ic < 
1£ B(i, Jd) .W. DISTANCE, 7= =) ther 
1£ Bll) .CURRENT sees” = (B(1, (J- 


abs (B(1I,J) .W. DISTANCE) ))) (GURRENS Rees 
+abs(B(1I,J).W.WEIGHT)) then 
B(1,J) .CURRENT COST -= (ba. 
abs (B(1I,J) .W.DISTANCE) ) ) .CURRENT_COST 
+abs(B(I,J) .W.WEIGHT) ) ; 
Bl. J) spree erien c= “WES 
Siqiel skies 
end sar - 
endy ii 


end FIND_MIN; 


—_— www wees aes ese 


begin 
while O@=. 9)! Bieo. 
MARK = = SEs 0ay 


GZ 


GOGUNTER <=. 1; 


[o@eiimee ion, NOP Loop 
EOm@mesthe t .. NOP “loop 
Peer Neal oreNCk s=1; 
bel) 2 DlSeoaniCk ==] 
Bie obi otaNncre <=] 
Big) Webbe lANCE. <= 1: 
BSAC alae: JU) he. 
end loop; 
end loop; 


Ooehee wr MebE —= IN FILE, NAME => "“ter.dat"); 
Evo ne, DIMENSTON OF MATRIX ee or 


Gb? 2 NOP) > 

PUT_LINE ("ENTER THE DELTA (the jumping value)."); 

EO eee | = ") Cee DELI) ; 

Puce oiNe VENTER THE VOLTA (the optimization tolerance) .") ; 
ve veto = ee Ger (VOL PA) 

Piao ieee Fite ke Pee sOURCE POINT !")-; 

Poe ene) Stoo @ee eur "SY = ")>-GET (SY); 

Put bikes ENTER The. COAL POINT !"); 

Omer en=e Ge) (GX) PUT( GY = ") -GET(GY) > 


-- CLOK function begins to compute the execution time. 
Soak eo Pe. = -ChOK ; 


cme mm ei e 


tow ROW in |. .NOP = loop 
Poco n 1 NOP gigas 
cet elvr ae (ROW, COl) CURRENT COsT) ; 
end loop; 
S90) lgele) em 
mee bOk LUNctVen finismes the computation of execution time. 
END Piva. = iChLOkK-: 


PU. CtOveteeriMe fO READ THE DATA FILE IS: peas 

PUT (END_TIME-START TIME, 4,4,0) ;new_line; 

START TIME = CLO; > 
--It determines the borders and calculates the weights of the 
edges... 
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vem ea ea i ei Ce Ce 


for LI in 2, -NeGeP eee 
for J in, | NOR rece 
CALSWEIGHE (iad). 
end loop; 
end loop; 


i ie ee ee 


--It makes the costs Max Mumber 1n Graders temice then: 


comparisons --for finding the memameneena. 
for I in l=. NOP oes 
for J in “Ll See NOP  eieoe 
B(lL,J).CURRENT COST 2== cc 0cGe 
B(l,J) OLD COST -:= *“BG@iyd) 2euLe et 2eo ie 
end loop; 


end loop; 


Ce rr Nd 


B(SX, SY) .CURRENTDZ¢COS)-- 0) fee 
loop 
A eaeea Dy Ee 
COZ = Ue 
1f A <= DELI then 
e@p 


exit when (A > DELI) ; 
if (J>=NOP) then exit; end 1f; 
1£((I=SX and J=SY) or (I=GX and J=GY)) and (MARK =TRUE) 


then 
MARK := FALSE;exit; end 1; 
MARK <=> 2RUE 
J =i oe ee: 
CO fue = he OC ea 
A -2= A + BL) -E. WELCH. 
Snide loop. 
Vi eOZe OVE hen 
B( lj J-2O2) -E-wErGHY =A = Bilsad) .B2WETen 
B(1ljJ-GO2Z) .E.DESTANGCE = =3 eer. 
B(1I,J) .W.WEIGHT -= Bip -COZi EE -Weuclie, 
BLL,U) .W. bist snes S= B(1,U-COZ) Be bist alee, 
end 1f; 
else 
Do 2S ad eee 
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ene ek - 
fon = NOP then 
[ee ee 5 
a Se he 
Simic: jf: 
exit when I = NOP + 1; 
ema loop; 


De i i i ee 


Ne 


Pie. obo lANCE = 1 then 
Zoe yd). DISTANCE + J; 


loop 
woe wp 
Sse Wem =. 2: 
Pied) ceo tae bh, 2s. =] 
Siglo wd wave LoreNnCe. <= —1; 
GCs = COSo+ 2; 
SiC eu Mere) oF 
else 
J. ®= Ajfi> 
end 1 £; 
Pe = NOP enen 
SEU te Ee aes 
Wis 
jelol sisae 


exit when I = NOP + 1; 
ena Loop; 


eee ewe eee 


emma ee eS 


A f=) bot wd) aN WE LGHT * 
GOz oO 
feel = se DBL Chen 
loop 
exit when (A > DELL) ; 
1f (I>=NOP) then exit; end 1 f; 
Pe hoevoreiced=cr or (l=Gx and J=GY)) and (MARK =TRUE) 
then 
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MARK := FALSE exit; encase 


MARK @s=" fave. 

aL = eae 

COZ == eOZ a. 

A s= A + B(T,Jd) -NeWElens 
end loop; 


Lf CO2)>. Ov ehen 

Bi T-€O2 3) NW ler 
B( i=CO2, J). NepistaNves 
Bil) s. Veer 
B(l,w) 2s -DistA ies 


A > Bly) NeWErGHT 
COM 

B( L-C0Z,7 0) ey Eerie 
B(L-€O77 od} ib hee ern, 


end if; 
else 
ali c=, 2 ae ir: 
end 1€£; 
1£ I = NOP then 
Ui <2 = ae 
de — a 
Veer wage % 
exit when J = NOP + 1; 
eric eer, 


ee 


1f Bi(1,J)  N-DIS TANCE =] 2 ehen 
4 += Bid) Nels ANGE oe 
iele) eo 
i= Jeele 
exit when I = Z; 
BGld aN DLs. 2Ner 
Bild) So. VEStANcE. = =a 
COS 3 = CO 
end loop; 
else 
| eee oll. 
end if; 
ie JE 
i = + 1; 
i — ale 
end if; 
exit when J = NOP + 1; 


iI 
I 
} 


NOP then 


— a 
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end loop; 


ee ee ee 


cme eee i 


OOS oS 0 3a — Ss ae a 
mon 1 ta lo... NOP leoep 
Maton cnceml=>oy) dnd mot (l=GxX and [=GY) then 
Ges: =) 0; 
Peek tl) Nee heraANe ky <1. Chen 
GO4) >=" CcoO4 = 1 
=iniclzalar 
iowa bf) ke PlolANer =< 1 then 
GO4- 4 = COAL ci: 
S101C ee 
feet ft) eo Pio TANCE < 1. then 
G@4 == CO4 = 1 
end if; 
cos il ware rANCe < 1 then 
CO4w = COdr ee, 1a 
end if; 
1f €CO4 >= 3 then 
Balieet nC hives PALSE : 
ieee) NeeiboraNnGe 7 =.) then 
BemieraesiS (heb) eNeDEiSTANCE)),1I)-S.DISTANCE := -1; 
Pai et DLSTANGE. «= -1% 
COs ==CeOse+ 1: 
end if; 
ieee ah ie anem (=. — 1 then 
poli ams ott) fs. Dlso LANCE) ))2Ww.DESTANCE := -1; 
Pie) woe oPANGh s= =) 
GOs P=COsa ee 
end if; 
it eeeh rs MP roTANGH /= —-1 then 
PG abso ola!) o. DLSTANCE))7J).N.DISTANCE ;:= -1; 
Pile Gwe wo olLANCE <= 2-1; 
GO3 ==COse.. 1. 
end if; 
Pee he Ww aDESrANGs /= —1 then 
Peele aicn Si i, lew. DISTANCE) )).B.DISTANCE := -—1; 
Eee ease ANG Re c=- =): 
COSTAE OS) sees ie 
end if; 


oy 


COS = GOs 


palol Ses 

end if; 

for J an i427)... (Ner Woes 
CO4 ean 


1f not(I=SxX and J=SY) and new 1l_@ ene ee 


1f Bl(l,d) .N,DISTANGE = etna 


GO4>:= Coa we. 

end if; 

1 B(1L,J) -E-DISTANCE = 1] then 
CO4 y= Ceca: 

end if; 

1f£ B(l,J) .S,DESTANCE = sl eeneu 
CO4 := C04 + 1; 

end if; 

1£ B(I,J) .W.DISTANCE < 1 then 
202 = OO a: 

end if; 

1£ CO4 >= 3 then 
joe Cllr ers OIE cchauece—, ars dLysych G 


1£° B(i,J) N. Bis PANGE Sy —)— eae 


B((I+abs(B(I,J) .N.DISTANCE) ) ,J) 
B( I,J )4AN-DIsTANnces == 2. 
CO3° 2 =COsr ae 

end if; 


lt Bil,J) PE PESTANCE 7—s eee 
Bll (J+abs (a(t, 2) eb eis rANcE ae 


B(il,J) -EB2DIsTANCE 2-2-1. 
CO3> *=Coum 
end if; 


1£ B(I,J).S. DISTANCE Weel then 


B( (T-abs(B{( LA). S-DisSTtaAnicl ear 
Bit.) Ge bis TANGE =.= 1. 
COs. =COs =k. 
end if; 
1£ B(l,J) .W.DISTANGE (= 5-1) ea1en 
Bil, (J=-abs(B(i A) Weiser ee 
Bl J) 2WebislAnce =... 5) 
GO34. =COnn ae 
end if; 
COs =" Cosy +e. 
end if; 
end af; 
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.W.DISTANCE 


W.DISTANCE 


.N.DISTANCE 


EB Ure taANnce 


then 


oa = 0 > 


Poenem—Sxnaqna I=sY) amd not(s=GX and I=GY) then 
ee hia es bisrANCe < 1 then 
Ge4= CO4 + 1; 
Siiceorr- 
Peete or OS PANCE = 1 then 
CO4e = co U4 s+ 1: 
end if; 
Pe Bid, ties. DlorANCe =< £ then 
COE — 1 eOr Bees pale 
end if; 
me 8), Ll) WebISTANCE = 21 then 
Col: = Coa, I- 
end if; 
Meme == oe einen 
Ea itn eie ye @e = FALSE: 
Peed) NePISPANCE /= =1 then 
Paw toabs (B40, 2) oN. DISTANGE)),1)'.S.DISTANCE := -1; 
Er ieee) oT ANGE wee 
COs v=COsrr 1 
end 1£:; 
Pees) Seen STAN@E )/ = —— iiehen 
Bidet itabs (bi, ie E.DMSTANCE))) .W.DISTANCE := -1; 
Bid, l) See it STANCE := =1; 
Goo -cOs = i: 
end i€f; 
Pewee ig nl) .5.DISTANG: ~—= —1 then 
Pwo doo on bho LANCE joie DLkSTANCE <= -—L- 
Pig il) So. otANGE =e — 1. 
GOs =COs 4 as 
end ie 
Rees) oP LS Al@e =f Ehen 
Silas te eae DP hSorANGe ey) ) ok DISTANCE := -L: 
EHee bs LANG B= il: 
GOS 22 CO3 + i: 
Gnade e- 
Coos. COs: +13 
end if; 
Sine 1b - 
end loop; 
end loop; 
Puce ees.) > 


——_ www see eee 


as: 


=-CoSt Minmimazat tenes 
while COUNTER > 0 loop 
SOQonyUae s= (bi 
FOr @ly nil NOE se ap 
for J in 12 NOP @ioce 
FIND MEIN ( Lone. 


VOLT := B(1,J) .OLD.COST —-Bi(1i,J) -CURES Tiaee. (- 
1f° VOLT —> VOmieAwenen 
COUNTER -:= COUNTER. 
else 
COUNTER += s@OUL TER: 
Sma ioee 
Bil,J) -OLD_COST <= Bi iy) -euUrR REN ee su 
end loop; 
end loop; 


end loop; 


SS 


loop 
PUSH(S,B(GX, GY) ,GURKENIZeC sae 
PUSH. (Sl B(GX, GY jie Pre rerron) 
SO’ -=B (GX GY) DIRECTO: 
if ~SQO =/"NORTH” then 
PUSH (S;,B(GX, GY )aNeDrEStANnGEs 
GX :=oGR+B (GA Gy) ere reer. 
elsif SQ = "EAST " then 
PUSH(S,B(GX, GY) .Bebisraner 
GY t= GY+B(GA,GY) .E. BPIsSTAaANce- 
elsif SQ) = 3 S@Uth Semen 
PUSH (S, B( GX, GY) 25 DistAliGrs 


GX := GX-B(GX,GY) .S.DISTANCE; 
elsif SQ = "WEST " then 
PUSH(S,B(GX, GY) .W. DISTANCE) ; 
GY <= GY-B(GxX, CGY Sy pier anGe- 
else 
ex1t; 
end if; 
exit when GX =SX and GY =SY ; 
end lo6ep; 
PUT iNew DISTANCE GOs DIRE @ REO hay .- 
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AUPE iy JUIN (a aie a ae 
loop 

PORIS Ewan h(E). PUT a 

Poo fe i BUTI E) PUT (" baie 


POS (Sh E) : 
ieee New then 
Poke cOUTH” ) -new line: 
elotresel = “BAST ”" then 
PUT ("WEST") ;new_line; 
ets igh eae] "SOUTH! St hen 
PUT ("NORTH") ;new_line; 
Sletrwei = §“"“WEST " then 
PUL YEAST” ) -mevw line- 
fae a8 © 


Scie ewoen S. LATEST =0; 
ania Loop: 
CLOSE (INE); 
--CLOK function finishes the computation of execution time. 
ENP TIME --= CLOK-; 
te TLaG TiMe TO EXECUTE THE PROGRAM IS : ae 
PUT (END_TIME-START_TIME,4,4,0);new_line; 


Ele "RUN ONE MORE TIME :y(es) or n(o) :") ;GET(Q); 
Meroe “on” then. exit; 

end if; 

end loop; 

end MAIN?2; 


al 


APPENDIX C: PARALLEL ALGORITHM SOURCE CODE 


--Title : PROJQ.ADA 
--Author —~ @ENGIZ sea) 
--Date 207067 92 
--Revised : 04/10/92 
-~=Course = FHESYsS 


-~-Compiler: ALSYS ADA 
--Description:Reads data from file,input starting and goal 
-- points,finds the minimum cost path. 
with TEXT_IO, COMMON, CHANNELS; 
use COMMON ; 
procedure PROJO 1s 
package INTEGER_INOUT is new TEXT_IO.INTEGER_IO(INTEGER) ; 
package FLOAT_INOUT is new TEXT_IO.FLOAT_IO(FLOAT) ; 
use FLOAT _INOUT, INTEGER. INGUT; 


INE : TEAT  TOSP LE Bs hy ee: 

B > GRID? GR: GRE pe rom 

counter : integer :=1; 

-- communication channels that are used 
OutToMars ; CHANNELS . CHANNEL REF == 


CHANNELS .OUTSPARAMERE RSs 2. )5 
InFromMars : CHANNELS .CHANNEL_REF 
CHANNELS .IN_PARAMETERS (2); 


mmm em em eee eee 


--This procedure calculates the weights on edges and builds 
up an hyphotetical 
--wall for processes to stay in the terrain... 


cme meme eee wee 


procedure CAL_WEIGHT (I,J : in INTEGER) is 


begin 
ic bee LG, then 
Bila) NawelGHt =e 
B({l;J).N.Bistance. |=) 
B(i+i,J7) . CURRENT IGCCr -—_o0CG; 
else 


B(l,J) .N.WEIGHT :=1+ ABS (Bil7d) CURE Ei i2Ces a 
Bil+l ie CURREIITE 2605 2 is 
end if; 
Tied = Logie eer 


WA 


Pewee rh WEIGHT := -1; 


el ee Pie raNCR s= -1;> 
Eee +i) CURRENT COST :=10000: 
else 
Pwd Peano c= LrARS{B( Lr we CURRENT COST - 
Peni). CURRENT COST) ; 
end if; 
iets = Ae Ae hen) 
Bat) oma BGHT <= -1; 
Bb ioeees  SLSMANCE == -1; 
Mole i AJ we UNMeiiieGeat =<=10000% 
else 


Peo Veer > — i ABSA Bi I,J).CURRENT_ COST - 
Eee) )  CURRENT_COST) ; 


enc 1k: 
Vee. = | then 

Bolte ee lenr ——s—1- 

El Soe bom ANC Kh = 1 

Bie g—-)) CURRENT s@0OST :=10000; 
else 


Eugenes). Wi WEEGHE :=1-eABS(Bil,J) .CURRENT_COST - 
B(I,J-1) .CURRENT_COST) ; 
end if; 
ena CAL WHIGHT- 
-- This part finds the minimum cost between current node and 
mes sf Our 
--neighbors (north, east, south,west) 
procedure FIND_MIN (I,J : in INTEGER) is 
begin 
sis CURRENT COST > (eal we CURRENT COST 
+abs(B(I,J).N.WEIGHT)) then 
Pere ne Onn tl ecogl b.wt) ) - CURRENT COST 
+abs(B(I,J).N.WEIGHT) ); 
NORTH” ; 


Bieta) ._DERBeCT ION 
end 1f; 
Pei eeURRENG scoot > Cie CURRENT eeGs ? 
+abs(B(I,J).E.WEIGHT)) then 
CS py CURRENT COST 
+abs(B(I,J).E.WEIGHT) ); 
Weed. is 


Pile eURe EN COST 


Ea j= DIREGTION 
end 1 f; 


v2 


1f Bil,J) (CURRENT @¢es 2s (Bt od) CURR Ree oar 
+abs(B(I,J).S.WEIGHT)) then 
B(1I,J).CURRENT_COST :=33 ea Seles mee = 
+abs(SB( Lo PeeeVibiGhtia:: 
Bro) Di REC imine “SOUTEs 
ened ant: 
1f BCU ld CURRENT Ie@Cs = (Bil ,J=!) (CURR E (ieee. 
+abs(B(I,J).W.WEIGHT)) then 


B(I,J) ,CURRENT_COS?T 32742 he |] Cah aren oan 
tabs (B(I,J) .W.WEIGHT) ) ; 
Bile pie eee Sa A Witol . 7s 
end if; 


end FIND_MIN; 
begin 

TEXT ITO.CPEN (INE ,MODE => DEX T20C2 ee hee 
"Cereal we) 4; 

TEXT 10O.PUT_LINE ("THE BDIMBNStON Cr Me tee =. 
INTEGERSINGUL- Gels. 

TEXT_IO.PUT_LINE ("THE NUMBER@Cr PROGEsSs CRS s=aen 
INTEGER INOGUTSGER 2. 

TEXT _IO.PUT_LINE (“ENTER THE OPTIMIZATION ICU ER (cian 

INTEGER_INOUT.GET (VOLTA) ; 

TEXT _IO.PUT_LINE ("ENTER DHE cCOURGE 2el it aa. 


TEXT_1LO.PUT_LINE( {3% =9)3)7, TNTEGER TNC UL 2 ae oe, 
TEXT_IO.PUT_LINE( SY =) See Teeen Ss NCU ere Ga 
TEXT IO. PUT_LINE | VENDTERS TE eee eo ere 
TEXT_IO.PUT_LINE ("GX = “) ; INTEGERS INOUT GEL ae 
TEXT_10.PUT_LINE ("GY = “) 7 INTEGERE NOU EIGER, 


meee eee 


— meme ees ss es eS es 


for ROW 1m. 1. 02 toop 
jong (TOIL) ahigh oil es, OP Lesye 
INTEGER_INOUT.{GEY (INP, BtROy, Cole CUReD ieee, 
end leep, 
end loop; 


—~ — ae ee ee eee eee eee ese esl 


www ww wi ees ee eS 


for 1 inwieeess cer 
for J lin ia Oeleor 
GRP =e ly 
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DATA_IO.WRITE (Out ToMars,GRI) ; 
end loop; 
end loop; 
--It determines the borders and calculates the weights of the 
edges.. 
mea t im 6:.104L00p 
poOmwoam Ls. k0- Loop 
CAL WEIGHT (I,J); 
end sloop; 
ence OOD? 
--It makes the costs max number in order to use them in 
Comparisons for finding 
ele oman MUM... 


mcm mw mw i 


E@ieek tin G.2h0 loop 


momen: Liweeeees (0 Loop 
Eel  CURRENPeCOstT = — 10000; 
Pee Ok DeaeOst ==) 8=b(l,0)scURRENT COST; 
eiotel - liayepe) 
ene loop: 


--This part sends dim of matrix,no of 

mis@eeesscors, VOlta,sourceand goal points. 

Peo Rae Nooo = 76k Er. OLD COST: =P;,-GRE.N.WEIGHT: =VOLTA; 
eee We boh tl aon GR Leo. WEILGHT:=SY>GRILCW.WEEGHT:=GA;GRI.E.DIS 
mANCE > =GY ; 

DATA_IO.WRITE (Out ToMars,GRI) ; 


meme rec mmm eee eee eee ee 


mmm meee as a ee 


Box, oY) CURRENT COST 
while COUNTER > 0 loop 
SouNreR = 0)- 
Pome die. .-0) Vboop 
IeOue qd) alin Ih see Ihe) a exeye 
heel. =VOmt nen 
Ghee = erhe a) 2 
DATA_IO.WRITE (Out ToMars,GRI) ; 
DATA_IO.READ(InFromMars,GRI)} ; 


I 
© 


fis 


B( i=l, 0). 2 =26Ee, 


end if; 
F ENDS Milian a = 
VOLT := B(1I,J).OLDLCOST —-Bi{ id) -eUFEe) | taeecu 
1£ VOLT = VO@—gAg@enen 
COUNTER .=.COUN TER 
enc “sere 
B(I,J) .OLD_COST <= Bild) CURRE ies. a. 
end loop; 
end <loep- 
integer_inout.. put (ceunmter) ; text s1toepue BO 0 
GRI.CURRENT [60ST = =Veouhia 


DATA..1O°. WRITE (OUCTOMars). Grin 
DATA_IO.READ(InFromMars, GRIT); 
COUNTER := GRE -CURREN i 2e et; 
integer Inoue pues  coumeen)- 
texb 10.put lane ("22324 2 2 oP 
end wleoe-; 
for fins 225. loop 
for J inviloe lOloece 
DATA_IO.READ( InbPromMars, "GR )- 
Jain Mie te SIL 
end. loop; 
end loop; 


loop 
PUSH(S, B( GK GY) CURRENT =GGsis 
PUSH] (S1,B(GxX7GY)— DERECTIGI ON 
5O -7=B(GxA GY) SDR Ee Lee 


1£ SO = "NORTH eaen 
PUSH(S;B(GX, GY) Sib tS TANCE 
GA se. GAT) 

elsif SOQ = "EAST " then 
PUSH (S73 (GA GY) 2 bop Po rAalie rs. 
GY eae 

elsif SQ = "SOUTH" then 
PUSH(S,B(GX, GY eo - Dl cTAnGE = 
GK 22 4GxX-L; 

elsift SO = ‘WES Semen 
PUSH (5, BGA, GY eee rs lal ere 
GY :=3Gy—i: 
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else 


exit; 

end if; 
Sole. WhenmGx =S5X end GY =SY -; 
=e Oop: 
mex) LO. PUT LINE (" DAE SYUANNeye Seowk Deore LON |) ~ 
iene TOG ENE (™ — ~=-== === ---- 8 --------- ee 
loop 


See ore LEGER. INOUL. PUT(E) »=TEXT.10O.PUT .(" ia 
POEs 2); INTEGER sIN@wWT.PUT(E);TEXT_IO;.PUT (" ble 


Berl Sl, ble 
feel = "NORTH" then 

heii OnPOiee SOUTH") ;TEXT TO.new line; 
Seti Racer  " then 

PEAT LOMEUT Ss 2WEST") ;TEXT_IO.new_line; 
Seni eh = "SOUTH" then 

Reet hoe Pur WUNORTE”)-TEXT.1O.new_line;: 
elsif El = "WEST " then 

Ca corO.rlT A bAStT )> TEXT 1LO.new line; 
end if; 
exit when S.LATEST =0; 

Sor Loop; 


ie LO CLOSE CINE ) ; 
Sma PROG; 


vy 


--Title : PROJ1.ADA 


--Author Te CENGE AZ EB rer 
--Date A One oe 
--~Revised : 04/10/92 

=~ eles elas alls 
-~-Compiler: MERIDIAN ADA 
--Description: 


with TEXT_IO, COMMON, CHANNELS; 

use COMMON; 

procedure PROJ1 is 

-- communication channels that are used 

OutToEarth : CHANNELS.CHANNEL_REF := CHANNELS.OUT_PARAMETERS 
(Zane 

InFromEarth : CHANNELS.CHANNEL_REF := CHANNELS.IN_PARAMETERS 
(203 

B : (GRID? GRI : GRID_POINT; 
counter,counterl - INtegeia am, 


cme ee i 


procedure CAL WEIGHT (I,J : in INTEGER) is 


begin 
Tf i. -= 10 then 
Bf) he WElee te —. 
Bl 1l,J)2NeDiS Paver: —sak 
B(14+1,J). CURRENT 2CGSE = =! 0000. 
else 


B(i,J) .N.WEIGHT s=l+ ABS(Bi(1i, 7) CURRE lee. 
Bit+l wa) sCURRENTZ¢6es i, 


end if; 
ae ie — IO Veiga. 
Bild) EWE lens... 
Bil, Jd) 2 blo wANeRa) — aoe, 
Bil ,d+ lL) ,@uURREND. Ces lucie 
else 
B(i,;J) .E.WEIGHT :;= I+ABS(B( 1a) )2eue. =) tees o a. 
B(il,J+)) VQuURREN Eee Ose, 
Srelcl shies 
fea een 
B(i;d) 5S. WEIGH. =: =s— 1. 
Bii,J).S -SrstAnes = ——— 
BCr-l, J) -CURBENT COC: = 00 ca: 
else 


B(1I,J).S.WEHELGHT :=l+ ABS(Bi( 1,7) -CURRERNI eos (eres 
1,J) CURRENIEGes 1); 
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end if; 
If J = 1 then 


Pele ewo WE LGHT <= =: 

Biwi Se ohANGCE 2:= =—1-> 

Pee See eCURRE NT seOST -=10000; 
else 


Peewee weet he: 1 ABS (B(T,J) .CURRENT COST - B(I,J- 
me CURRENT COST); 
eneeif; 
end CAL WEIGHT; 
procedure FIND_MIN (I,J : in INTEGER) is 
begin 
Pe tO CURRENT COST > (ee ed) .CURRENT VeOST 
+abs(B(1I,J) .N.WEIGHT)) then 
aie rN heel w=" (Bi i+ 7) CURRENT COST 
+abs (B(I,J) .N.WEIGHT) ) ; 
ORR 


lo IC, .tJ) | eB Seev@ ILI) | 
end 1f; 
fe Bit Ji. CURRENT COST > Cea osCURRENT COST 
+abs(B(I,J).E.WEIGHT)) then 
(eee CURRENT COST 
Faos(S( Ind) Se Weer) ); 
al aiy ord Yair 


Bay eae ORnReNT [COST 


Bilao PREC TLON 
end if; 
1f B(I,J) .CURRENT_COST > Vee) | CURRENT COST 
+abs(B(I,J).S.WEIGHT)) then 
Diver ceo ie Ooi. = UP hi) CURRENT COST 
+abs(B(1I,J) .S.WEIGHT) ); 
YOOULEs 


Bold) . DIREGTION 
end if; 
ol nw CURRENT COST > (Sil, ) > CURRENT. GoOsT 
+abs(B(I,J) .W.WEIGHT)) then 
(Sie) See CURRENT COST 
+abs (B(I,J) .W.WEIGHT) ) ; 
WES bee’ 


Es CURRENT COST 


lev 10, OLE OAR ICON) 
end 1f; 
end FIND_MIN; 
begin 
Soule iil, .5 Loop 
EG ih ich Ny Ome ereye, 
DATA_IO.READ(InFromEarth, GRI) ; 


7s 


Billed) SSG r ie 
ence lecre; 
end loop; 
--It determines the borders and calculates the weights of the 
edges... 
for 1 an 1. 4 2a seeue@er 
FOr Win, he sO Seep 
CAL WEIGHT(I,J); 
end loop; 
end Weep: 
--It makes the costs max number in order to use them in 
comparisons for finding 
== the mMinrmumeens 
for I in 1 ...57lcep 
fOr J dt eeleen. FOMleces 
Bil,J).CURRENT COST. =e Ren ior 
B(1l,J).ObD-CGS2 -= Belg) euURE a wee sun 
end loop; 
end loop; 
--This part sends dim of matrix,no of 
proccessors,volta,sourceand goal points. 
DATA_IO.READ(InFromEarth,GRI); 
N:=GRI .CURRENT COST; P:=GRI.OLD COST; VOLU2.* =6Gr eae 
SX:=GRI.E.WEIGHT; SY :=GR®. S.WEILGHT ; Gx: =GRE_W-WeiGhl- GY. -Gr 
,DESTANGE. 


mmm wwe eee ees ses Ss 


B(SX, SY) JCURRENT= COs. — c- 
while COUNTER > 0 loop 
COUNRER=—  0- 
FOr oi sins. 1s. 7S loa 
for. in 1.4 OS iloep 
Li iss 3s) Steen 
DATA_IO.READ (InFromEarth,GR1I) ; 
Brel gd) 2 =" GE. 
GRE s=98 (150) ) + 
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Pera eto- vWRETe(OucToRarth,GRI) ; 
eng: : 

Peeve mrorn (1 4d); 

eee ee ie OLD ICOST —-B(1, J) seURRENT. COST; 

Pe Ol Men OLA “then 

SCoOOnwaRs: =. COUNTER +1; 

Smicuait - 

ie de Onmmee. terror Bit) -CURRENT COST; 
end loop; 

end loop; 

DATA IO.READ(InFromEBarth, GRI) ; 

@eounwer ll) 2= GRE?CURRENT COST; 

Pm OumceG= 0) and (Gounterl=0) then 
Ge eURREMmeCOs. == 0; 
DATA_IO.WRITE(OutToEarth, GRI); 

else 
COUNTER (a=se- 

GRI.CURRENT_COST :=1; 
Proto Vet Vem TORArEh, GRE): 
end if; 
ema Loop 
mote 2n L..5 loop 

Meme) Im | .710 loop 

Owe = B11) ) > 

PATASLO.WRITE (Out ToBarth,GRI) ; 

end loop; 

end loop; 
end PROJ1; 
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# File: makefile 

= ‘make help" to print) eperon (ee 

tt 

= Complete development cycle: 

# make family -- makes Ada family and library 
directories 

= make -- compiles, links, configures source 
# make run -= run bootable code 

MODE == 

PROC y=ase 


OPTS-= 7S(MODE) /Go Poe 


# make the executable code 
main.btl: mainh.c$(PROC)S (MODE) proejin cs (PROC) > Mone, 
main.pgm 
@ echo EXPECT 1 WARNING... 
1conf /s main.pgm 
@ TES Vit ie oedel 


mainh.cs$ (PROC)S (MODE): proj0.o pre UM Eo ekOeC > tere 
merger.t$(PROC)$(MODE) mainh.t$ (PROC) $ (MODE) 
pga iy cs (Me mmigisiatiaes Jo cle 


proj)0.02 common ada crejeeaca 
ada invoke proj0.inv,yes 


proj0h.t${PROC)S (MODE). proi0h2.tax prejUhvoce 
eECGaMm oi OFS) ~PLojUbeoee 


Dro7Un2 tax: preqUnz eee 
occam /ta /x, pre7 Un. oce 


merger.t$(PROC)$ (MODE): merger.occ 
occam $(OPTS) merger.occ 


mainh.t$(PROC)S$S (MODE): mainh.occ 
occam S (OPTS) marnheoce 


projlh.csS (PROC) $ (MODE): proj lo pred mie ot Oe ts Cnr 


1link projlh.t$(PROC)$(MODE) projl.o adarts8.1lib 
NOSti0. lias voeceanice ls 
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projl.o: common.ada projl.ada 
Doc eKe pro] l..inv,yes 


Peewee s (PROC) S (MODE) : projlh2.tax projlh.occ 
Seca OFS) projlh.occ 


pm@yqlne tax: proj lh2.occ 
Seca. (fax DrO7ln2 occ 


= 
# MISC. 
2 
help: 
@ echo Make arguments: 
@ echo make - make from top level down 
@ echo make -n [opt] - display but don't execute 
commands 
@ echo make *.0 - make Ada object 
@ echo make help ~ display this list 
@echo make clean - delete all files except source 
@ echo make run - run bootable program 
@ echo make check —ONeCCwiGransoputen topOLOgy 
@ echo make family - make Ada family and library 
directories 
clean: 
eel * 278? 
Glel* bax 
dele 
gel * asc 
del. ~.btl 
del test_lib\adalib.* 
rad test_lib 
del test_fam\adafam.* 
rd test_fam 
Te DQ E: 
1server /sb main.btl 
eneck: 
check /r 
family: 


ada invoke family.inv,yes 
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-- File: main.pgm 

#INCLUDE "hostio.inc" 
#INCLUDE “lamkacdeie sae. 
PROTOCOL PASS IS INR |[5)ByYTe 


#USE "mainh.c8s" 
fOSE: voword hecacy 


CHAN OF PASS Mars2Earth, Earth2Mars: 
CHAN “OF SP FromFiler. Ton mien. 


PLACED PAR 
PROCESSOR 0 T8& 

PLACE FromPiler "AT “biak02 in 

PLACE ToFiler AP AnkKO out: 

PLACE MarsZ2Farth AT link2.in: 

PLACE Barth2Mars AT link2.out: 

[13225000] IN vs l- 

main.harness (FromFiler, ToFiler, MarsZEarth, BartnZMarse 
wslL) 


PROCESSOR aks 


PLACE Barth2Mars AT Pinker: 
PLACE Mars2Earth AT linkO.out: 


[1230000] Nm wee: 
projl.harness (Mars2Earth, Earth2Mars, ws2) 
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-- File main.lnk 

-- Purpose: File list for ilink 

mainh.t8s 

merger.t8s 

mestic. lib 

occam8s.1lib 

MmeBen Omer os DlLOyU ,omee@arecomlib hostio.lib occam8s.lib } 


eri les fami ly.inv 
family.new test_fam, overwrite=yes 
lib(family=test_fam) .new test_lib,overwrite=yes 


eis: oro) 0. inv 

default.compile library=test_lib 

complile common.ada 

compile proj0.ada 

default.bind library=test_lib, level=bind,warning=no 
Dineen O]., Obj ect—"pro ]).c" ,entry point="pro)]0.program" 


niles Droge. 1ny 

default .compile library=test_lib 

Semolle proil.ada 

default .bind library=test_lib, level=bind, warning=no 
meA@more) bremece="Dreyl.o",entry point— prog il .program" 


oi) 


-- File: mainh.occ 
#OPTION "AGNVW" 
#EINGLUDE. “hes: oausmen 
PROTOCOL PASS 2S) Wi (se ee 
PROC main.harness (CHAN OF SP Frome iater Sten. 
CHAN OF PASS Mars2Earth, Earth2Mars, 
[J INT FreeMemory) 


#USE “hestio. lib: 


FUSE oe DcOmoMl se os. 
#USE "merger.té8s" 


[ARG HAN OF ANY Debug: 

[2]CHAN OF SP FromAda, ToAda: 

CHAN OF BOOL StopDebug, StopMultiplexor: 
SES 


PAR 


-- A multiplexor to combine the debug and normal output. 
so.multiplexor (FromFiler, ToFiler, FromAda, ToAda, 
StopMultiplexor) 


-- A debug channel merger. 
debug.merger (ToAda[0], FromAda[0], Debug, StopDebug) 


-- A process to invoke the sleve program. 
ws IS FreeMemory: 
SEQ 
proj0.harness (FromAda[1i], ToAda[1], Debug[0], 
Mars2Earth, Earth2Mars, ws) 
StopDebug ! FALSE 
StopMultiplexor ! FALSE 


so.exit (FromFiler, ToFiler, sps.success) 
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-- File: merger.occ 


FOPTION “AGNVW" 
PC wUes “hostao.a nc" 


PROC debug.merger (CHAN OF SP FromFiler, 
Lica AN Or ANY Debug, 
CHAN OF BOOL Stop) 


Hoo oHnOSstio.laib" 


-- A debug channel merger and blocker. 


Veuemax.caeoug 1S 20: 
Voienumoer.Oneaebug 1S SIZE Debug: 


PiGeline 71ndex: 

ise BYTE line. butter: 

BYTE values, rc: 

BOOL running, reset, Ss: 
{max.debug]BOOL mask: 

VAG BYTE laine.feed 1S 10 (BYTE) : 


Sy a @) 

SEOnt = U0 POR mumber Of .debug 
mask{i] := TRUE 

aUhaUG LOG) 8 BS IONS 

reset := FALSE 

Paanieom unde. = © 

WHILE running 
PRI ALT 


ALT 1 = 0 FOR number.of.debug 
mask[i] & Debug[i] ? value 
0) 
alee ts 
value = line.feed 
SEO 


TORseo 


-- Send the complete line. 


Seepucs (Frombi ler, 


Horuler, spice. stdout, 


[line. buffer FROM 0 FOR line.index], r) 


line.index := 0 

mask [1] := FALSE 

heset == TRUE 
TRUE 
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SEQ 
-- Add character tee! ine 


line.buffer[line.index] := value 
line.index := line vandex; a! 
reset & SKIP 
Sn© 
reset := "Ramo: 
SEQ 1 = 0 FOR Htiteer Jer cenug 
Ties Halen: =e ele, 
SE ODm Es 
running = PARSE 
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oie ore), Un.occ 


#OPTION "“AGNVW" 

MPNGGUDE. “nost lo. i1c * 

Prevoecon PASS IS INT; [5] BYTE 

PROC proj0.harness (CHAN OF SP FromAda, ToAda, 
CHAN OF ANY Debug, 
CHAN OF PASS Mars2Earth, Earth2Mars, 
[J INT FreeMemory ) 


#IMPORT "projOh2.tax" 


[1] INT dummy.ws: 
wsl IS FreeMemory: 
Po eNT <neprogram: 
fo | ENT eue. pregram: 
Son ®, 
-- Set up vector of pointers to channels. 
tmorseqmemeel: -= MOSTNEG: INT -- not used 
LOAD.INPUT.CHANNEL (in.program[1], ToAda) 
Loree itneraANNEh (ineeregram{2], Mars2Earth) 
LOAD.OUTPUT.CHANNEL (out.program[0], Debug) 
LOAD .OUTPUT.CHANNEL (out.program[1], FromAda) 
HeAD oCUlTPU TT VOHANNEE (Out .euegram(2|, Earth2Mars) 
-- Invoke the Ada program. 
-- Assumes the entry point name has been changed to 
ieee) 0. orogram™. 
pum PGe@selmei wel, iM orogram, Out .program, Gummy.ws) 


== rile: PROJOH2Z.ccc 
ROErION "ARV" 


BeOe prey” program {{])INT wsl, in, out, ws2) 
Peo O OjeINT a: 
SEC 
SKIP 


89 


-- File: projlh.occe 


#OPTION "AGNVW" 
+#INCLUDE: “hostie sia]. 
PROTOCOL MASS Ss iL > exe 


PROC projl.harness (CHAN OF PASS Mars2Earth, Earth2Mars, 
[] INT FreeMemory) 


#IMPORT “promlnzaeca.. 


LD) INT dummy ews: 
wsl IS FreeMemory: 
(3) ENT “in = program. 
[3 ENDL? Oue. ere Gieam- 


SEQ 
-~- Set up vector of pointers to channels. 
in,program( 0] V = MOST: Gare -- not used 
in.programil|] >= MOSTNEGsaG -- standard i/o not used 
LOAD. INPUT.CHANNEL (in.program[2], Earth2Mars) 
out .program[0] := MOSTNEG INT -- standard i/o not used 
out.program[1] := MOSTNEG INT -- standard 1/o not used 


LOAD .OUTPUT.CHANNEL (Gut. pregranmi eliaseZsameu, 
-- Invoke the Ada program. 
-- Assumes the entry point name has been changed to 


"Dro erecranue 
projl.program (wsl, in.program, Sut -pregtam, sami ee 


=— Pale? prog ln scee 


#OPTION "AEV" 


PROC projl<program ( [| IND Wei euler. sc?) 
[22010000 |i ae 
SEQ 
SKIP 
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